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SUMMARY 


It  is  well  known  that  vehicles  traveling  at  high  speed  excite  the  air 
to  high  temperature,  resulting  in  dissociation  and  ionization,  so  that 
the  air  properties  deviate  considerably  from  those  of  a  simple  gas.  Some 
effects  of  these  reactions  are  considered  in  this  Report  in  relation  to 
the  thermodynamic  properties  of  the  air,  the  transport  coefficients,  and 
the  diffusion  of  heat. 

Methods  are  presented  for  calculating  the  thermodynamic  properties  of 
air  with  good  accuracy  up  to  15,000°  K  and  for  pressures  from  10"*  to 
10^  atmospheres. 


SOMMAIRE 


Ainsi  qu’il  est  blen  connu,  le  vol  a  grande  vitesse  de  tout  engln 
provoque  des  temperatures  elevens  de  I’air,  donnant  lieu  en  consequence 
aux  phenomfenes  de  la  dissociation  et  de  1’ ionisation,  de  sorte  que  les 
caracteristlques  de  I’air  different  sensiblement  de  celles  d’ un  gaz  simple. 
II  s*  agit  dans  ce  rapport  de  certains  effets  resultant  de  ces  reactions 
conslderes  en  fonction  des  propriet4s  thermodynamiques  de  I’air,  des 
coefficients  de  transport  et  de  la  diffusion  de  chaleur. 

L’  auteur  pr^sente  des  precedes  de  calcul  des  caracteristlques  thermo¬ 
dynamiques  de  1’  air,  avec  une  bonne  precision  pour  des  temperatures 
s’ elevant  a  15.000°K  et  des  pressions  comprises  entre  10‘*  et  10^  atm. 
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NOTATION 


R  universal  gas  constant 

p  pressure 

p  density 

T  temperature 

molecular  weight  at  standard  conditions 
Z  total  number  of  moles  from  one  mol  of  initially  undlssoclated  air 

fraction  of  oxygen  molecules  dissociated 
fraction  of  nitrogen  molecules  dissociated 
fraction  of  atoms  ionized 

Kj  equilibrium  constant  for  reaction  0^  20 

equilibrium  constant  for  reaction  ^  2N 

population  -  weighted  average  of  equilibrium  constants  for  reactions 
0  0**^  Y  e“  and  N  +  e“ 

mol  fraction  of  chemical  species 

H  enthalpy  per  mol 

S  entropy  per  mol 

partition  function  for  i^^  species 

y  ratio  of  specific  heats 

c  speed  of  sound 

V  velocity  of  flow 

velocity  of  flow  at  initial  conditions 

X  distance 

t  time 

r  distance  between  atoms 

cr  collision  diameter 


V 


mass  diffusion  diameter 


►  constants 

j 

electron  charge 
polarlzabil ity 
total  charge  density 
viscosity 

mean  molecular  velocity  of  i^*^  species 
mean  free  path  of  i^*^  species 
particle  concentration 

partial  coefficient  of  thermal  conductivity 
heat  flux  potential 
thermal  diffusivity  (=  k/pC^) 
coefficient  of  heat  capacity 


Prandtl  number 


THERMODYNAMIC  AND  TRANSPORT  PROPERTIES  OP 
HIGH-TEMPERATURE  AIR 

C,  Frederick  Hansen* 


1.  INTRODUCTION 

It  is  axiomatic  that  the  science  of  aerodynamics  must  be  based  on  a  good  understand¬ 
ing  of  the  medium  through  ^vhich  vehicles  fly.  At  subsonic  and  low  supersonic  speeds 
near  the  earth's  surface,  it  is  generally  sufficient  to  assume  that  the  atmosphere 
behaves  as  an  ideal  gas  and  that  the  ratio  of  specific  heats  is  a  constant  with  a  value 
about  1,4,  However,  it  is  well  known  that  vehicles  traveling  at  high  speed  excite  the 
air  to  high  temperature  with  the  result  that  air  properties  deviate  considerably  from 
those  of  a  simple  gas.  For  example^  Figure  1  shows  the  major  chemical  reactions  which 
are  produced  in  the  stagnation  regions  of  vehicles  as  a  function  of  altitude  and  velo¬ 
city.  At  about  3000  ft/sec,  the  vibrational  energy  of  air  molecules  is  excited. 

Oxygen  dissociation  begins  near  7000  ft/sec  and  nitrogen  dissociation  occurs  at  velo¬ 
cities  greater  than  15,000  ft/sec.  Single  ionization  of  atoms  becomes  dominant  near 
escape  velocity.  The  dissociation  reactions  are  encouraged  by  low  pressure,  of  course, 
and  the  altitude  dependence  shown  in  Figure  1  is  a  consequence  of  this  fact.  Some 
effects  of  these  chemical  reactions  will  be  considered  pertaining  to  the  thermodynamic 
properties  of  air,  to  some  of  the  transport  coefficients,  and  to  the  diffusion  of  heat 
through  gases. 

2.  EQUILIBRIUM  THERMODYNAMIC  PROPERTIES  OF  AIR 

The  equilibrium  thermodynamic  properties  of  gases  can  be  calculated  with  good  con¬ 
fidence,  provided  the  energy  levels  of  the  gas  particles  are  accurately  known,  and  a 
number  of  such  calculations  have  been  tabulated  for  air.  Perhaps  the  most  complete 
and  exacting  calculations  are  those  performed  by  Gilmore^  and  by  Hilsenrath  and 
Beckett^,  Both  of  these  calculations  are  very  exact  in  the  sense  that  they  take  into 
account  many  of  the  minor  chemical  components  of  air  and  small  corrections  to  the 
partition  functions  of  the  molecules.  Logan  and  Treanor^  have  also  prepared  tables 
of  air  properties  for  which  all  of  the  important  components  are  considered.  These 
precise  calculations  involve  multiple  iteration  procedures  such  that  the  labor  involved 
would  be  discouraging  indeed,  if  it  were  not  for  electronic  computers. 

The  content  of  the  tabulations  in  References  2  and  3  has  been  summed  up  graphically 
in  the  Mol  Her  diagrams  prepared  by  Feldman^  and  by  Moeckel^.  Using  such  graphs  one 
can  solve  a  number  of  important  aerodynamic  problems,  such  as  finding  the  properties 
of  isentroplc  expansions  or  of  non-isentropic  shock  discontinuities  (see,  e.g. ,  Ref.  4, 
6,  7,  or  8). 

For  a  number  of  reasons,  it  is  desirable  to  have  approximate,  analytic  solutions 
for  the  thermodynamic  properties  of  air  in  addition  to  the  precise  solutions  mentioned 
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above.  From  an  inspection  of  the  tabulated  equilibrium  concentrations  of  the  chemical 
species  in  air\  it  is  immediately  apparent  that  there  are  but  four  reactions  of  major 
importance  over  a  wide  range  of  pressure  and  temperature.  These  are  the  dissociation 
of  molecular  oxygen  and  of  molecular  nitrogen,  and  the  ionization  of  atomic  oxygen  and 
atomic  nitrogen:-  , 


0^ 

-  20 

(1) 

-  2N 

(2) 

0 

^  o'"  +  e* 

(3) 

N 

--  N'^  +  e' 

(4) 

With  one  exception,  all  other  reactions  yield  concentrations  which  are  of  the  order  of 
0.1%,  or  less,  of  the  major  components  given  by  these  reactions.  The  exception  is  the 
formation  of  nitric  oxide,  which  may  become  a  sizable  minor  component  of  air  at  high 
pressures  (see  Ref,  1).  However,  the  heat  of  formation  of  NO  is  very  nearly  the 
average  of  that  for  0^  and  ,  so  the  production  of  NO  does  not  appreciably 
change  the  balance  between  atoms  and  molecules.  As  a  result,  the  thermodynamic  func¬ 
tions  for  air  are  not  strongly  influenced  by  this  reaction. 

Some  distinctive  features  of  the  foregoing  reactions  are  convenientiy  summarized  on 
a  graph  of  the  compressibility  as  a  function  of  temperature  (Fig,  2),  The  compressibi 
lity  is  defined  by  the  gas  law 


_P 

P 


(5) 


where  R  is  the  universal  gas  constant,  p,  p  and  T  are,  respectively,  the  pressure, 
density  and  temperature  of  the  gas,  and  is  the  molecular  weight  at  standard  con¬ 

ditions.  To  the  approximation  that  the  particles  obey  the  ideal  gas  law,  Z  represents 
the  total  number  of  moles  formed  from  a  mol  of  initially  undissociated  air.  Since  air 
contains  about  20%  oxygen,  the  compressibility  approaches  1.2  when  oxygen  dissociation 
is  complete.  It  increases  further  to  about  2.0  when  nitrogen  dissociation  finishes  the 
conversion  of  molecules  into  atoms.  Single  ionization  of  the  atoms  doubles  the  number 
of  gas  particles  again,  so  that  Z  approaches  4.0  when  these  reactions  are  complete. 
Note  that  the  slope  of  Z  gets  very  small  between  reactions,  indicating  that  one 
reaction  is  essentially  complete  before  the  next  begins.  This  follows  from  the  fact 
that  the  dissociation  energy  of  nitrogen  (9,76  ev)  is  nearly  twice  the  dissociation 
energy  of  oxygen  (5.08  ev),  while  the  ionization  energy  of  nitrogen  atoms  (14.54  ev) 
and  oxygen  atoms  (13.61  ev)  is  nearly  three  times  as  great.  Thus  the  different  reac¬ 
tions  proceed  in  different  intervals  of  temperature,  and  for  purposes  of  approximation 
it  can  be  assumed  that  they  are  independent  of  one  another.  The  two  ionization  reac¬ 
tions  occur  simultaneously  but  with  nearly  the  same  energy  changes,  so  for  this  case 
it  is  assumed  that,  once  the  air  becomes  dissociated,  all  atoms  constitute  a  single 
species  which  has  the  population-weighted  average  properties  of  the  nitrogen  and 
oxygen  atoms. 


We  shall  review  briefly  some  of  the  essential  relations  involved  in  calculating 
the  approximate  thermodynamic  properties  of  air.  Further  details  may  be  found  in 
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Reference  9.  If  is  the  fraction  of  oxygen  molecules  which  are  dissociated, 

the  fraction  of  nitrogen  molecules  dissociated,  and  the  fraction  of  atoms  which 
are  ionized,  then  the  compressibility  is 

Z  =  1  +  +  2^3  (6) 

The  reactions  are  assumed  to  be  completely  independent  and,  to  the  order  of  approxima* 
tion  being  considered,  it  seems  adequate  to  take  the  ratio  of  nitrogen  to  oxygen  in 
air  as  4  to  1.  Then  one  obtains 


/  1 

-0.8  +  yo.64  +  0.8(l  +  — 


2(1 


-0.4  +  / 0.16  +  3.84  1  + 


f  4p 
2  1  +  — 
V  K,  . 


4p 


1 

V/(1  +  P/Kg) 


(7a) 


(7b) 


(7c) 


where  and  are,  respectively,  the  equilibrium  constants  for  reactions  (1)  and 

(2)  in  terms  of  partial  pressures;  K3  is  the  population-weighted  average  of  the 
equilibrium  constants  for  reactions  (3)  and  (4).  The  equilibrium  constants  are  calcu¬ 
lated  from  partition  functions  which  include  only  those  energy  levels  up  to  6  kT  (8.76 
ev  at  ISpOOO^K).  The  energy,  entropy  and  specific  heats  of  the  components  are  also 
calculated  from  these  abbreviated  partition  functions. 

The  form  of  Equations  (7a)  through  (7c)  is  suitable  for  constant  pressure  calcula¬ 
tions.  For  constant  density  calculations,  these  equations  become 


-1  +  / 1  + 


3.2pRT 


8/ljRT 
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At  this  point  it  may  be  remarked  that  the  formation  of  NO  can  be  included  in  the 
closed  form  approximations  if  one  is  willing  to  go  to  the  trouble  of  solving  cubic 
equations  rather  than  the  foregoing  quadratic  expressions. 


The  mol  fractions  x.  of  the  chemical  species  are  given  by 


0.2  -  e. 


'2  Z 

0.8  -  e. 


2  Z 

2^  *  0.463 
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26^  -  1.663 
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(9a) 


(9b) 


(9c) 


(9d) 


(9e) 


and  the  enthalpy  and  entropy  for  the  mixture  are  found  by  a  simple  summation  over  all 
the  species 


H 

RT 


d(ln  0^) 
6 (In  T) 


(10) 
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d(ln  Q^) 
6  (In  T) 


x^dn  x^) 


(11) 


The  symbols  H  and  S  refer  to  the  enthalpy  and  entropy  per  mol  of  particles  and 
0^  Is  the  partition  function  for  the  species  at  the  standard  state  pressure  p^. 

Note  that  the  and  O'**  fractions  should  be  considered  separately  for  the  term 
x^(ln  x^)  which  appear  in  the  entropy  (Eq.  11),  though  they  may  be  combined  elsewhere. 


The  accuracy  of  the  foregoing  approximations  is  indicated  in  Figures  3(a)  through 
3(d).  Figure  3(a)  shows  the  percentage  discrepancy  between  the  compressibility  given 
by  Equation  (6)  and  the  solutions  of  Hilsenratb  and  Beckett^  for  densities  of  10^.  10*^ 
and  10*^  Amagat,  Figures  3(b)  and  3(c)  show  the  same  comparison  for  enthalpy  and  entropy. 
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It  can  be  seen  that  the  discrepancies  are  generally  less  than  5%,  sufficient  for  many 
engineering  purposes.  The  largest  discrepancies  occur  at  high  densities  where  there 
is  some  overlapping  of  the  reactions  and  where  the  formation  of  nitric  oxide  should  be 
taken  into  account  for  greater  accuracy. 


The  specific  heats  of  air  may  be  found  by  taking  the  derivatives  of  energy  and 
enthalpy.  The  approximations  lose  accuracy  with  each  successive  order  derivative,  of 
course,  but  still  the  agreement  is  generally  within  10%  of  the  iteration  solutions 
performed  by  Logan^,  Fortunately,  one  does  not  usually  use  the  specific  heats  directly 
in  aerodynamic  calcnlatinns  an^nvay,  but  rather  the  ratio  of  specific  heats,  y.  As 
shown  in  Figure  3(d),  the  approximations  preserve  the  accuracy  of  this  quantity  within 
a  few  percent. 

The  speed  of  sound,  c,  will  be  defined  by 


c 


2 


(12) 


which  may  be  transformed  to 
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Using  this  quantity,  one  can  evaluate  the  integral 


(13) 


I 


(14) 


which  is  useful  for  determining  the  Rieraann  invariants  v  ±  I  for  one*dlmensional , 
isentropic  flow.  The  change  in  velocity  along  the  characteristic  directions  in  such 
gas  flow  is  given^^  by 


V  -  -  ±  ~  1)  (15) 

where  the  subscript  i  refers  to  initial  conditions  and  the  positive  sign  corresponds 
to  the  positive  characteristic  direction 

dx 

—  -v+c  (16a) 

dt 

and  the  negative  sign  to  the  negative  characteristic  direction 


dx 

dt 


V  -  c 


(16b) 
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in  the  time-distance  (x,t)  plane.  The  integration  of  Equation  (14)  is  along  a 
constant  isentrope,  and  the  solutions  are  performed  by  iteration  with  an  electronic 
computer  inasmuch  as  entropy  is  not  one  of  the  independent  variables  in  the  analytic 
expressions  for  the  gas  properties.  Where  very  small  changes  in  flow  properties  are 
desired,  as  in  the  construction  of  a  fine  net  characteristics  solution,  it  is  more 
convenient  to  work  with  just  the  integrand  of  Equation  (14).  Small  changes  in  flow 
velocity  are  approximately  given  by 


V 


p\  Ap 


cp/  p 


(17) 


where  the  algebraic  sign  of  the  pressure  change  Ap  corresponds  to  the  sign  of  the 
characteristic  direction.  The  quantity 


I 


is  plotted  in  Figure  4  normalized  by  the  ideal  air  value 


I* 


(18a) 


(18b) 


where  7  is  taken  equal  to  1.4.  The  ratios  are  given  as  functions  of  temperature  for 
isentropes  ZS/R  equal  to  30,  40,  50,  60,  '70,  80,  90  and  100.  Values  for  intermediate 
entropies  may  be  estimated  reasonably  well  by  interpolation.  Numerical  values  of  I/I* 
and  I  are  tabulated  in  Reference  11. 


3.  TRANSPORT  PROPERTIES  OF  AIR 

Next  consider  some  additional  properties  of  air  needed  for  aerodynamic  calculations, 
namely  the  transport  properties  viscosity  and  thermal  conductivity.  In  estimating 
these  properties,  one  is  faced  with  two  alternatives.  The  transport  coefficients  may 
be  calculated  by  the  fairly  rigorous  but  complex  method  of  Chapman  and  Enskog^^  or  by 
the  simple  but  approximate  kinetic  theory  of  hard,  elastic  spheres^^.  If  good  esti¬ 
mates  of  the  interparticle  potentials  were  available,  so  that  the  Enskog  collision 
Integrals  could  be  determined  accurately,  the  choice  would  immediately  fall  on  the 
former  method,  of  course.  However,  at  present  the  potentials  are  uncertain  enough  so 
that  it  is  doubtful  whether  the  inherent  accuracy  of  the  more  rigorous  method  could 
be  realized.  Inasmuch  as  encouraging  progress  Is  being  made  on  quantum  mechanical 
calculations  of  the  interparticle  potentials  (see,  c.g. ,  Refs,  14  through  17),  it  has 
been  decided  to  defer  the  labor  of  making  the  rigorous  solutions  until  all  of  the 
important  potentials  are  available.  In  the  meantime,  some  estimates  based  on  the 
simple  kinetic  theory  have  been  prepared.  It  may  be  noted  that  this  elementary 
approach  need  not  be  discounted  unduly,  for  it  does  give  reasonably  good  answere  where 
experimental  evidence  is  available  for  comparison. 

The  first  step  is  the  calculation  of  mean  free  paths  from  the  effective  collision 
diameters,  Hirschfelder  and  Eliason'®  have  found  that  for  a  wide  variety  of  different 
attractive  potentials,  the  collision  integral  method  corresponds  to  an  effective  hard 


7 


sphere  diameter  equal  lo  the  distance  where  the  Interparticle  potential  is  about  •0. 6 
kT,  In  the  case  of  repulsive  potentials,  the  effective  diameter  for  momentum  and 
energy  transport  is  the  distance  between  the  atoms  when  the  potential  is  about  0.9  kT* 
while  for  mass  diffusion  it  occurs  at  about  1.6  kT.  These  simple  criteria  might  be 
applied  if  only  the  potentials  were  known.  Unfortunately,  in  the  case  of  atom-atom 
and  atom- ion  collisions,  the  picture  is  complicated  by  the  fact  that  the  particles  may 
approach  one  another  along  any  one  of  a  number  of  different  potentials,  depending  on 
how  the  electron  spin  vectors  happen  to  add  up.  For  example.  Figure  5  shows  the 
qualitative  form  of  the  potentials  for  collisions  between  two  neutral,  ground-state 
nitrogen  atoms  as  a  function  of  the  distance  r  between  atoms^^.  The  lowest  lying 
potential,  ^2,  has  the  lowest  total  electron  spin  and  it  is  known  quantitatively  from 
the  vibrational  energy  levels  observed  spectroscopically  for  the  stable  ground  state 
molecule.  Since  the  higher  lying  potentials  are  not  known  quantitatively,  we  are 
forced  to  estimate  the  average  effective  size  from  the  single  known  potential.  For 
this  purpose,  it  has  been  assumed  that  the  average  effective  collision  diameter  is  the 
distance  between  the  atoms  when  U(^2)  equals  -kT  for  the  case  of  momentum  and  energy 
flux,  and  -2  kT  for  mass  diffusion.  This  procedure  undoubtedly  overestimates  the 
size  of  shallow  intermediate  potentials  (such  as  ^2  and  ^2  -  Fig.  5),  On  the 
other  hand  it  underestimates  the  relatively  long-range  forces  between  particles  in 
excited  states.  The  very  steep  repulsive  potentials  (such  as  ^2  of  Fig.  5)  are 
roughly  symmetrical  with  the  lowest  attractive  potential  over  the  long-range  part  of 
the  internuclear  separation  distance,  so  this  potential  is  probably  correctly  accounted 
for  to  a  first  approximation. 

To  be  more  specific,  the  following  formulas  were  used  for  calculating  collision 
cross  sections  (a  Sutherland  type  cross  section  is  used  for  oxygen  and  nitrogen 
molecules): 


(19) 


where  the  constant  C 
temperature,  is  3.14  x 
may  be  approximated  by 


is  112®  K  and  the  molecular  cross  section  at  very  high 
10“^^  cm^.  The  lowest  lying  atom-atom  and  atom-ion  potentials 
the  Morse  function 


U 

D 


-;3{r/r„-i)1^ 
1  -  e  ® 


1 


(20) 


whence  the  collision  diameter  for  momentum  and  energy  transfer,  a,  is  given  by 


a  1  /  /  kT 

-  =  I  .  -  In  ll  -  I  -  — 
P  \  i  D  , 


(21a) 


and  the  diameter  for  mass  diffusion,  cr' ,  is 


—  =  1  -  -  In  I  I 


2kT\ 

'•t) 


(21b) 


The  constants  r^,  JS,  and  D/k  for  the  various  types  of  collisions  are  taken  from 
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Herzberg^^  and  these  are  listed  in  the  following  table: 


Co  1 liding 
partic  les 

D/k 

(°K) 

a 

'■e 

(Angstroms ) 

0-0 

59,  000 

3.24 

1,207 

N-N 

113,200 

2.96 

1.094 

N-0 

75,400 

o  to 

O  4  X 

1.  151 

0-0'^ 

75.20C 

3.18 

1.123 

N-N* 

101,200 

2.94 

1.116 

For  electron-atom  collisions,  the  atom  becomes  polarized  by  the  approach  of  the 
electron  and  a  charge-dipole  type  of  interaction  results.  The  classical  interaction 
energy  of  the  induced  dipole  with  the  electron’s  field  is 

e^a 

U  = - -  (22) 

2r‘* 

where  e  is  the  electron  charge  and  a  is  the  polarizability.  The  polarizability  is 
taken  to  be  13.2  x  and  10.3  x  for  oxygen  and  nitrogen  atoms,  respect¬ 

ively.  The  quantum  mechanical  cross  sections  calculated  by  Hammerling,  Shine  and 
Kivel^*^  show  that  the  Hartree  electrostatic  field  and  the  exchange  integrals  both 
contribute  to  the  Interaction  energj'  as  well  as  the  polarization.  However,  they  find 
that  these  two  effects  roughly  compensate  one  another,  and  so  these  corrections  will 
be  omitted  for  the  approximation  being  considered  here. 

Finally,  the  ion-ion  and  ion-electron  cross  sections  are  calculated  from  the  simple 
Coulomb  potential  between  two  charges.  The  correction  for  the  charge  shielding  effect 
given  by  Cohen,  Spitzer  and  Routly^*^  makes  the  cross  sections  weakly  dependent  on 
density  as  follows: 


S  =  —  ]  In 


kT 


1  /2 


2e'^  \7m^ 


(23) 


where  n^  is  the  total  charge  density,  both  positive  and  negative. 

The  viscosity  was  calculated  from  the  simple  summation  formula  for  a  mixture  of 
hard  spherical  molecules-^: 


577-  V 

55  Zfi'-ih 


(24) 
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where  is  the  mean  molecular  velocity  of  the  species  and  the  mean  free  path 

is  given  by 


i 


(25) 


where  n  is  the  particle  concentration.  Equation  (25)  differs  from  Kennard' s  expres¬ 
sion  for  the  mean  free  path^^  to  account  for  the  fact  that  the  momentum  of  particle  i 
relative  to  particle  j  is  decreased  by  the  fraction  2(1  +  M^/Mj)"^  on  each  elastic 
collision  with  particle  j.  This  factor  multiplies  each  term  of  the  summation  in  the 
denominator  of  Kennard’ s  formula,  and  the  result  is  the  same  as  obtained  in  the  more 
rigorous  matrix  formulation  for  the  viscosity  of  mixtures^ ^  if  the  off-diagonal  ele¬ 
ments  of  the  matrix  are  set  equal  to  zero.  This  is  valid  in  a  first  approximation  at 
least,  as  the  off-diagonal  elements  are  small  compared  to  the  trace  elements  of  the 
matrix.  Buddenberg  and  Wllke^^  have  shown  that  experimental  data  on  the  viscosity  of 
mixtures  can  be  reproduced  by  empirically  adjusting  the  numerical  factor  2  before 
those  terms  of  the  summation  in  Equation  25  which  Involve  unlike  particles  (i  f  j). 
However,  in  view  of  the  other  approximations  Involved  in  the  present  analysis,  it  did 
not  seem  warranted  to  use  a  more  sophisticated  formulation  for  the  viscosity  of  the 
mixture  than  given  above.  The  estimated  viscosities  are  shown  as  a  function  of  tem¬ 
perature  for  pressures  of  10°.  10"^  and  10"“  atmospheres  in  Figure  6.  The  ordinate  is 
the  ratio  of  the  coefficient  given  by  Equation  (24)  to  the  value  given  by  the  Suther¬ 
land  formula: 


^0 


1,462  X  10"® 


112 

1  +  - 

T 


gro 

cm  sec 


(26) 


where  T  is  in  degrees  Kelvin.  Presumably,  Equation  (26)  is  about  the  viscosity  which 
the  gas  would  have  if  air  molecules  remained  inert.  The  effect  of  the  dissociation 
reactions  is  to  increase  slightly  the  ratio  m/Mq.  However,  the  charged  particles  have 
very  large  cross  sections  when  they  collide  with  one  another,  so  the  viscosity  is 
reduced  to  comparatively  low  values  at  the  higher  temperatures  where  ionization  becomes 
predominant. 


A  theory  for  the  thermal  conductivity  of  a  chemically  reacting  gas  was  first  out¬ 
lined  by  Nernst^®,  The  theory  is  used  here  in  the  form  subsequently  developed  by 
Hirschfelder^**,  in  which  the  energy  transfer  through  the  gas  is  treated  in  two  inde¬ 
pendent  parts.  The  first  mode  of  energy  transfer  is  by  molecular  collisions  just  as 
for  ordinary  non-reactive  gases.  The  second  mode  of  energy  transfer  is  by  diffusion 
of  the  molecular  species  and  the  reactions  which  occur  as  the  gas  tends  to  maintain 
itself  in  chemical  equilibrium.  For  the  first  mode,  use  is  made  of  Eucken’ s  assump¬ 
tion^®  that  the  Internal  energy  is  distributed  among  the  gas  particles  independently 
of  the  velocity  distribution,  so  that  the  partial  coefficient  of  thermal  conductivity 
takes  the  form 
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Hirschfelder  formulates  the  reactive  contribution  to  thermal  conductivity  in  terms  of 
the  multicomponent  diffusion  coefficients  which  are  somewhat  difficult  to  estimate. 
However,  Butler  and  Brokaw^^  have  shown  that  the  reactive  coefficient  of  conductivity 
may  be  expressed  in  terms  of  the  ordinary  binary  diffusion  coefficients, 


^reactive 


nk 


fd(ln  K) 
\d(ln  T), 


X  .D 


i'^ij 


ajXi) 


(28) 


where  the  are  the  stoichiometric  coefficients  of  the  components  when  the 

chemical  reaction  is  written  in  the  form 


=  0 


(29) 


The  diffusion  coefficients  are  calculated  by  elementary  kinetic  theory  as  follows^^: 


J_  ^  sA  _^V 
Djj  3  \^T 


(30) 


The  total  coefficient  of  thermal  conductivity  given  by  Hirschfelder^"*  is  just  the  sum 
of  the  inert  and  reactive  parts: 


^  ^inert  ^  ^reactive 

The  estimated  coefficients  of  thermal  conductivity  are  shown  for  pressures  of  10^, 
10"^  and  atmospheres  as  a  function  of  temperature  in  Figure  7.  Again  the 

ordinate  is  normalized  by  an  inert  gas  value  given  by 


19  R 


(32) 


The  thermal  conductivity  goes  through  pronounced  maxima  at  temperatures  where  the 
chemical  activity  is  intense.  As  shown  in  Figure  7,  the  coefficient  may  become  one 
to  two  orders  of  magnitude  lar^-pr  the  inert  gas  value.  The  maxima  occur  at  about 
the  same  temperature  where  the  compressibility  changes  most  rapidly  with  temperature 
(Fig.  2)  and  therefore  where  the  specific  heats  are  a  maximum.  In  fact,  Hirschfelder^** 
has  noted  that  the  thermal  conductivity  should  be  approximately  proportional  to  the 
specific  heat  for  the  chemically  reacting  gas  just  as  it  is  for  the  inert  gas  (Eq.  27). 
It  follows  that  the  Prandtl  number  should  not  be  grossly  affected  by  the  reactions,  at 
least  up  to  the  temperature  where  the  ionizatlcn  reactions  occur.  The  highly  ionized 
gas  is  a  mixture  of  very  light-weight  particles  (electrons)  with  heavy  particles  (ions) 
and  the  effect  of  this  mass  difference  is  to  increase  greatly  the  ratio  of  thermal 
conductivity  to  viscosity.  Thus,  if  the  ionized  gas  is  everywhere  in  thermal,  equil i- 
brium,  it  should  have  a  very  small  Prandtl  number,  of  the  order  of  10"^  (see  Ref.  9). 
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Some  very  approximate  assumptions  are  used  in  the  foregoing  calculations  of  the 
transport  properties,  and  it  is  obviously  desirable  to  check  them  with  experiment  in 
some  manner.  It  seems  difficult  to  obtain  the  transport  coefficients  directly  by 
experiment  at  the  very  high  temperatures  involved,  but  it  has  been  possible  to  evaluate 
a  quantity  called  the  heat  flux  potential  by  combining  some  results  of  shock  tube 
experiments  with  the  theory  of  heat  diffusion  through  gases.  The  heat  flux  potential 
is  defined  as 


T 

^  ^  dT  (33) 

Jq 

so  it  represents  a  transport  property.  This  integral  has  the  quality  of  a  potential, 
since  the  heat  flux  at  any  point  in  a  medium  is  just  grad  <^,  The  heat  flux  potential 
is  the  natural  parameter  to  use  as  the  dependent  variable  in  problems  involving  heat 
diffusion  through  media  where  the  thermal  conductivity  Is  a  function  of  temperature, 
as  it  is  for  gases.  In  terms  of  0,  the  heat  diffusion  equation  takes  the  familiar  . 
form 


—  .aV^0  =  0  (34) 

ot 

where  a  is  the  thermal  diffusivlty,  Note  that  the  dlffusivity  cannot  be 

considered  a  constant  as  is  usually  done  in  the  classical  solutions  of  the  heat  dif¬ 
fusion  equation  (see  Ref.  27,  for  example). 

Consider  now  the  one-dimensional  heat  flow  through  a  semi-infinite  gas  medium.  We 
shall  restrict  ourselves  to  the  case  where  the  initial  conditions  and  the  boundary 
conditions  are  constants;  that  is,  the  initial  potential  throughout  the  gas  is  the 
constant  designated  by  0^^,  and  the  boundary  value  is  the  constant  0^  (see  Fig.  8). 

The  diffusivities  are  and  a^,  respectively.  Normalized  coordinates  will  be  used 

such  that  and  are  both  taken  to  be  unity.  The  results  can  easily  be  general¬ 

ized,  of  course,  to  account  for  arbitrary  units  of  these  boundary  conditions. 

The  idealization  of  constant  initial  and  boundary  conditions  will  be  approximately 
realized  in  a  physical  situation  where  the  reservoir  is  a  slab  of  material  with  a  very 
large  heat  capacity  and  a  large  thermal  conductivity,  such  as  a  metal.  Then  the  wall 
can  soak  up  heat  fast  enough  to  maintain  nearly  constant  conditions  at  the  interface. 
The  slab  might  be  suddenly  immersed  in  a  constant  temperature  gas,  for  example;  or, 
alternatively,  the  gas  might  be  heated  suddenly  by  a  plane  shock  wave  reflecting  from 
the  solid  interface. 

Boltzmann^®  has  shown,  with  perfect  generality,  that  the  solution  to  this  problem 
may  be  expressed  as  a  function  of  x//t.  Therefore,  the  time  and  distance  variables 
are  transformed  to  the  single  dimensionless  parameter 


X 

■/(4a^t) 


(35) 


where  is  the  diffusivlty  of  the  gas  at  the  boundary;  is  unity,  of  course, 

in  the  present  normalized  coordinate  system,  but  this  is  not  an  essential  feature  of 
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the  transformation.  In  terms  of  this  parameter,  y,  the  partial  differential  equation 
of  heat  conduction  becomes  the  dimensionless,  total  differential  equation 

d^(^  64> 

a  — r  +  2y  —  =  0  (36) 

dy^  dy 

The  factor  4  in  the  parameter  y  (Eq.  35)  is  an  arbitrary  stretching  factor.  It 
is  chosen  merely  so  that  the  solutions  to  Equation  (36)  reduce  to  the  usual  error 
function  form  when  the  diffuslvlty  is  constant,  that  is,  when  a  equals  The 

merit  in  using  the  normalized  form  for  the  potential  (that  Is,  <p/(p^)  is  that  one 
set  of  integrations  for  Equation  (36)  will  suffice  for  all  possible  boundary  conditions. 
All  that  is  required  before  one  proceeds  with  the  integration  is  that  the  diffusivity, 
a,  be  evaluated.  Two  cases  will  be  considered,  one  in  which  the  gas  is  ideal  and 
inert,  and  the  second  in  which  the  gas  is  in  local  equilibrium  but  is  chemically  active 
with  a  large  heat  of  reaction. 

According  to  the  kinetic  theory  of  inert  gases^^,  the  coefficient  of  thermal  con¬ 
ductivity  is  approximately  proportional  to  a  power  b  of  the  temperature,  where  b 
is  close  to  Vz.  Thus  the  integral  of  thermal  conductivity  is  proportional  to  the 
(b+l)^*’  power  of  temperature: 

cp  -  J K  a  (37) 

Now  the  thermal  diffusivity  equals  the  conductivity  divided  by  the  heat  capacity 
per  unit  volume,  CyO.  Since  the  density  p  is  inversely  proportional  to  temperature, 
for  an  ideal  gas,  the  diffusivity  is  proportional  to  the  (b+1)^^  power  of  temperature 
also; 


a 


•pb+i 


(38) 


It  follows  that  diffusivity  is  proportional  to  the  potential  and  in  the  present 
normalized  coordinates  this  means  that 


a  -  0  (39) 

Sometimes  a  linear  relation  between  a  and  4>  with  a  finite  Intercept  will  best 
fit  measured  values,  but  this  can  always  be  transformed  to  the  direct  proportionality 
of  Equation  (30)  by  appropriate  adjustment  of  the  lower  limit  of  the  integral, 
Equation  (33),  Generally,  this  limit  will  be  close  to  absolute  zero. 


It  may  be  remarked  that  even  if  one  accounts  for  the  variations  of  heat  capacity 
which  occur  in  a  real  gas,  the  preceding  relations  are  valid  to  the  order  of  approxi¬ 
mation  that  Prandtl  number  is  a  constant.  This  can  be  seen  from  the  fact  that  diffu¬ 
sivity  is  just  the  kinematic  viscosity  divided  by  the  Prandtl  number: 


_  p/p 

a  -  - 

Pr 


(40) 


and  the  kinematic  viscosity  is  again  approximately  proportional  to  the  (b+1)^^  power  of 
temperature,  at  least  in  gases  composed  of  neutral  particles. 
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When  the  result  of  Equation  (39)  is  used,  the  heat  diffusion  equation  takes  on  the 
inocuous  looking  form 

6^d)  di 

4^—^2y—  =  0  (41) 

dy'^  dy 

Analytic  solutions  to  Equation  (41)  are  known,  but  none  which  satisfies  the  boundary 
conditions,  namely: 

<^(0)  =  1  (12a) 

lim  -  <4, 

r-oc  (42b) 

This  is  somewhat  frustrating  in  view  of  the  deceptively  simple  appearance  of  the  non¬ 
linear  Equation  (41).  However,  it  is  relatively  easy  to  integrate  the  equation 
numerically,  starting  from  a  given  value  of  the  boundary  derivative,  (^'  =  (dt^/dy)^, 

and  terminating  as  the  solution  asymptotically  approaches  a  limit,  The  solutions 

are  something  like  error  functions  stretched  slightly  out  of  shape.  Three  of  these 
solutions  are  shown  in  Figure  9  as  a  function  of  The  solutions  ^  are  shown 

divided  by  the  error  functions  1  +  (0^  -  l)erf  y/Vcp^.  The  0  functions  rise  more 
steeply  than  the  error  functions  near  the  orgin,  go  through  a  maximum  deviation,  then 
approach  the  error  function  as  a  limit. 

The  value  of  <4o  is  uniquely  related  to  the  derivative  at  the  origin,  4>q-  This 
relation  is  shown  in  Figure  10.  The  derivative  is  equivalent  to  a  dimensionless 

heat  flux  at  the  boundary,  while  the  normalizing  function,  1  +  0',  is  the  value 

which  4^^  would  have  if  diffusivity  were  a  constant.  We  shall  return  to  this  relation 
between  and  4>q  after  considering  the  case  of  the  chemically  reacting  gas. 

A  chemical  reaction  in  the  gas  behaves  like  a  reservoir  which  soaks  up  heat  as 
temperature  is  increased,  and  liberates  heat  when  temperature  drops.  Consequently, 
the  specific  heat  is  very  large  if  the  heat  of  reaction  is  large  compared  to  RT. 
Hlrschfelder^**  has  shown  that  the  coefficient  of  thermal  conductivity  for  gases  in 
local  equilibrium  is  also  very  large;  in  fact,  it  is  approximately  proportional  to 
the  specific  heat.  Thus,  the  Prandtl  number  is  relatively  constant  and  the  diffusi¬ 
vity  is  a  function  of  temperature  which  is  not  greatly  affected  by  the  chemical 
reaction.  The  Integral  of  thermal  conductivity,  on  the  other  hand,  is  greatly  Increased 
as  a  result  of  the  reaction.  This  situation  is  illustrated  in  Figure  11. 

The  diffusivity  increases  linearly  with  4  up  to  the  point  where  the  chemical 
reaction  occurs,  there  it  flattens  out  until  the  reaction  is  about  complete,  then  it 
increases  again  in  a  more  or  less  linear  manner.  The  solid  lines  are  the  previously 
discussed  theoretical  estimates  for  air  in  which  oxygen  dissociation  occurs.  The 
pressure  dependence  of  the  curve  for  4>  Isss  than  24^  has  been  removed  by  normal¬ 
izing  both  ordinate  and  abscissa  with  the  factor 


4 


(43) 


b+i 
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where  is  the  boundary  temperature,  and  is  the  temperature  where  the  thermal 

conductivity  Is  a  maximum,  that  is,  where  the  diffuslvity  is  most  nearly  Independent 
of  the  conductivity  integral.  The  temperature  may  be  calculated  from  Equation  (28). 

It  should  be  noted  that  the  calculation  of  the  relation  between  diffusivlty  and  the 
conductivity  integral  does  not  depend  on  precise  numerical  values  of  the  transport 
coefficients.  For  example,  although  the  estimate  shown  for  air  in  Figure  11  is  based 
on  simple  kinetic  theory  and  very  approximate  collision  cross  sections,  the  corrections 
introduced  by  more  exact  calculations  affect  both  the  diffusivlty  and  the  conductivity 
simultaneously  in  such  a  way  that  the  functional  relation  between  the  two  Is  maintained. 


The  relation  for  reacting  gases  in  general  will  be  similar  to  that  for  air.  The 
dashed  curve  on  Figure  11  is  a  limit  which  is  approached  as  the  heat  of  reaction  becomes 
very  large.  This  limiting  relation  has  been  used  in  the  integrations  which  follow,  that 
is,  a  is  assumed  to  equal  up  to  the  point  and  thereafter  is  taken  to  be 

constant.  Up  to  the  point  then,  the  solutions  are  the  stretched*out  error  func¬ 
tions  which  were  discussed  earlier  and  which  are  shown  in  Figure  9.  At  0^,  these 
solutions  are  Joined  by  the  solution  for  constant  diffusivlty  which  has  a  matching 
slope  at  the  junction^  (dc/^/dy)^: 


4> 


exp  ~ 


(44) 


In  Figure  12,  the  limit  of  the  potential  is  shown  for  a  reacting  gas  as  a 

function  of  the  boundary  heat  flux,  The  solution  now  depends  on  the  value  of 

4i^,  which  in  turn  depends  on  the  boundary  temperature,  the  pressure,  and  the  chemical 
reaction  being  considered.  For  example,  in  the  case  of  air  and  a  boundary  temperature 
of  0°  C,  the  values  20,  40  and  60  for  correspond  to  pressures  of  about  lO""*,  10° 

and  10^  atmospheres,  respectively.  The  solutions  for  values  of  corresponding  to 
other  pressures  and  different  boundary  temperatures  can  be  obtained  fairly  accurately 
by  interpolation  between  the  numerically  integrated  solutions  such  as  the  curves 
plotted  in  Figure  12.  The  inert  gas  solution  (Fig.  10)  corresponds  to  4>^  equal  to 
infinity  (Pig.  12).  It  should  be  noted  that  these  solutions  can  only  be  used  up  to 
the  point  where  the  chemical  reaction  goes  to  completion.  Beyond  this,  the  increase 
in  diffusivlty  must  again  be  taken  into  account.  The  domain  of  validity  is  a  charac¬ 
teristic  of  each  specific  reaction.  H'ithln  this  limitation,  the  curves  of  Figure  12 
apply  to  chemically  reacting  gases  in  general. 

The  above  solutions  have  been  put  to  use,  in  conjunction  with  experiment,  to  evalu¬ 
ate  the  heat  flux  potential  as  a  function  of  temperature^’ '^°.  In  this  experiment,  a 
plane  shock  wave  is  reflected  from  the  end  wall  of  a  shocic  tube.  The  heat  flux  lu  Llie 
wall  is  measured,  and  this  fixes  the  abscissa  for  the  graphical  solution  on  Figure  12, 
The  temperature  at  the  wall  and  the  pressure  of  the  gas  after  the  shock  reflection 
determine  the  value  of  (see  Eq.  43),  Hence  the  value  of  the  ordinate  in  Figure 

12  is  fixed,  and  the  value  of  4:^  is  determined.  The  temperature  of  the  gas  associa¬ 
ted  with  this  heat  flux  potential  <4,  is  taken  to  be  the  equilibrium  temperature  after 
the  shock  reflection,  and  this  is  a  known  function  of  the  measured  shock  velocity. 

This  procedure,  of  course,  implies  that  the  experiments  are  conducted  under  conditions 
where  the  chemical  relaxation  times  are  short  compared  to  the  test  interval. 
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Results  of  the  experiments  and  calculations  are  shown  in  Figure  13.  The  heat  flux 
potential  is  shown  divided  by  the  inert  gas  value,  where  the  coefficient  b  is  taken 
to  be  1/2 


^inert 


2  . 

^  ~  ^  3.2  X  10  ^ 

3  »  ‘\tJ 


cal 

cm  sec 


(45) 


and  where  T  is  in  degrees  Kelvin.  Thus  the  ratio  unity  if  the 

coefficient  of  thermal  conductivity  is  proportional  to  the  l/2-power  of  temperature. 
This  appears  lu  be  approximately  true  up  to  the  temperature  where  oxygen  dissociation 
begins.  At  the  pressures  involved  in  the  present  shock  tube  experiments,  this  occurred 
at  about  2500^K.  At  higher  temperatures,  the  heat  flux  potential  becomes  more  than 
twice  as  large  as  the  inert  gas  value  as  a  result  of  the  dissociation  of  oxygen.  The 
solid  curve  on  Figure  13  shows  the  theoretical  heat  flux  potential,  based  on  the  esti¬ 
mates  for  the  coefficients  of  thermal  conductivity.  The  agreement  between  the  experi¬ 
ment  and  theory  is  about  as  good  as  the  scatter  in  the  data  and  is  somewhat  better  than 
the  uncertainty  expected  of  the  theory.  The  small  deviation  from  inert  gas  behavior 
which  occurs  just  before  the  dissociation  reaction  has  been  explained^^  as  the  result 
of  small  increases  in  the  thermal  conductivity  integral  due  to  vibrational  energy 
excitation  and  to  formation  and  ionization  of  nitric  oxide.  The  theory  and  experiment 
appear  to  agree  on  the  large  change  due  to  oxygen  dissociation,  and  it  is  believed  that 
the  theoretical  estimates  may  represent  a  reasonable  approximation  to  use  for  some 
engineering  purposes.  However,  it  is  clear  that  the  experimental  results  give  merely 
a  necessary  but  not  sufficient  check  on  the  validity  of  the  approximate  methods  used 
to  evaluate  the  transport  coefficients.  It  may  be  noted  that  the  potential  shown  in 
Figure  13  is  not  a  constant  pressure  function  but  follows  the  particular  progression 
of  pressures  obtained  with  the  various  strength  shock  waves  produced  for  the  experi¬ 
ment.  The  theoretical  potentials  are  shown  in  Figure  14  for  temperatures  up  to  SOOO^K 
and  for  constant  pressures  from  10^  to  10*^  atmospheres.  Interpolations  can  be  made 
approximately  in  proportion  to  the  logarithm  of  pressure. 


4.  CONCLUDING  REMARKS 

To  summarize,  some  methods  have  been  outlined  by  which  thermodynamic  properties  of 
air  can  be  approximated  in  closed  form  with  an  accuracy  generally  somewhat  better  than 
5%  ^nr  temperatures  up  to  15,000‘^K  and  for  pressures  from  10*^  to  10^  atmospheres. 

The  solutions  are  analytic  in  terms  of  temperature  and  pressure  or  of  temperature  and 
density  as  the  independent  variables.  A  single  series  of  iterations  suffices  for 
calculation  of  constant  entropy  properties. 

The  equilibrium  mol  fractions  have  been  used  for  estimates  of  the  viscosity  and 
thermal  conductivity  for  air.  These  estimates  are  based  on  the  elementary  kinetic 
theory  of  colliding  elastic  spheres  where  the  size  of  the  sphere  is  allowed  to  vary 
with  temperature  so  as  to  account  for  some  functional  dependence  of  collision  cross 
section  on  velocity.  The  coefficient  of  thermal  conductivity  integrated  with  respect 
to  temperature  is  a  heat  flux  potential  and  is  found  to  be  the  natural  parameter  to 
use  as  the  dependent  variable  in  problems  involving  heat  diffusion  through  gases. 

This  potential  has  been  evaluated  by  experiment  for  air  up  to  5000°K,  including  the 
effect  of  oxygen  dissociation,  and  the  results  are  found  to  agree  reasonably  well  with 
the  theoretical  estimates. 
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Pig, 3  Deviations  in  the  approximate  solutions  for  thermodynamic  properties  of  air 
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Fig, 6  Viscosity  for  air 
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Fig. 12  Chemical  equilibrium  heat  diffusion 


Pig.  13 


Experimental  heat  flux  potential  for  air 


Pig. 14  Theoretical  heat  fJux  potential  for  air 


DISTRIBUTION 


Copies  of  AGARD  publications  may  be  obtained  in  the 
various  countries  at  the  addresses  given  below. 

On  peut  se  procurer  des  exeraplaires  des  publications 
de  1' AGARD  aux  adresses  suivantes. 


BELGIUM 

BEbGIQUE 


CANADA 


DENMARK 

DANEMARK 


FRANCE 


GERMANY 

ALLEMAGNE 


GREECE 

GRECE 


Centre  National  d' Etudes  et  de  Recherches 
A4ronautiques 

11,  rue  d'Egmont,  Bruxelles 

Director  of  Scientific  Information  Service 
Defense  Research  Board 
Department  of  National  Defense 
‘A’  Building,  Ottawa,  Ontario 

Military  Research  Board 
Defense  Staff 
Kastellet,  Copenhagen  0 

O.N.E.R.A.  (Direction) 

25,  Avenue  de  la  Division  Leclerc 
ChS,tillon-sous-Bagneux  (Seine) 

Wissenschaftliche  Gesellschait  fur  Luftfahrt 
Zentralstelle  der  Luftfahrtdokumeutation 
Munchen  64,  Flughafen 
Attn;  Dr.  H.  J.  Rautenberg 

Greek  National  Defense  General  Staff 

B.  MEO 

Athens 


ICELAND  Director  of  Aviation 

ISLANDS  c/o  Flugrad 

Reykj  avik 


ITALY 

ITALIE 


LUXEMBURG 

LUXEMBOURG 


Ufficio  del  Generale  Ispettore 
del  Genio  Aeronautico 
Ministero  Difesa  Aeronautica 
Roma 

Obtainable  through  Belgium 


NETHERLANDS  Netherlands  Delegation  to  AGARD 

PAYS  BAS  Michiel  de  Ruyterweg  10 

Delft 


NORWAY 

NORVEGE 


PORTUGAL 


TURKEY 

TURQUIE 


UNITED  KINGDOM 
ROYAUME  UNI 


UNITED  STATES 
ETATS  UNIS 


Mr.  0.  Blichner 

Norwegian  Defence  Research  Establishment 
Kjeller  per  Lillestrom 

Col,  J,A.  de  Almeida  Viama 
(Delegado  Naclonal  do  ‘AGARD’ ) 

Direc^ao  do  Service  de  Material  da  F.  A. 

Rua  da  Escola  Politecnica,  42 
Lisboa 

■Ministry  of  National  Defence 
Ankara 

Attn.  AGARD  National  Delegate 

Ministry  of  Aviation 
T. I. L. ,  Room  009A 
First  Avenue  House 
High  Hoi  born 
London  ff.  C,  1 

National  Aeronautics  and  Space  Administration 
(NASA) 

1520  H  Street,  N.W. 

Washington  25,  D.  C. 


Printed  by  Technical  Editing  and  Reprodvetion  Ltd 
95  Great  Portland  St.  London,  k\l. 


es  O)  I 

iH  W. 

aJ  fli  O 

O  G 

CO  m 

CO  O  ^  C!j 

x:  4J 

y  o  * 

TO  ^  q>  jj 

^  -C  « 

O  c  -iJ  0; 

^  O 


<y 

c:  -c  « 

4J  ;3 

nU  'H  'l-C 
Ui  O  ■■^ 
C  T3 
Cl  CO  z,^ 

£  .2  ^ 


r:  T3‘ 
x:  i  c 
+j  c*  ed 
o 


JL, 

cd  05 

1 

■i-i 

-rH 

Li 

Lc 

l. 

<d 

<iH  Cd 

0; 

O 

o 

X 

G 

05 

to 

<0 

X 

tl>  c 

c 

•G 

K  O 

<15 

Ol 

+-> 

o 

X 

+-> 

It 

+J 

(d 

o 

X 

cd 

-M 

o 

•M 

X 

Cd 

O  ” 

c 

X) 

0) 

o 

o 

X 

to 

<1-1  0) 

-r^ 

to 

-!-» 

<i-l 

* 

cd 

•.-< 

o 

j3  a> 

i  o  .5  ■ 

5  « 

^  U  ^ 

S  ^ 

c "  a 

<J  c  f  <y 
<D  a: 

:  cj 

<y 

,  e  c 

'  ■'^  o  ^ 


-!->  ■4->  y 

G  C  3  5 

O  ci  cd  z 
G  «  c  6 
a>  -M  o  X 

DS  <  t-  0 
0)  2 
O  x:  <  * 
Di  ^  Oi 

o  o  O  S 

<1:  z:  ^ 


(-  iti 

□.  c>  a 


U  •  <7i  i£> 
H  O  •-I  M 


M  c  •-' 

S  C  a 

to  ^  ii-i 

C  "o 

O  o  o 

O  E  o 


U  .  O  CO 

E-.  O  <M 


“  ^  T3 

Q>  <d  5 
to  2 

+J  (t 
*-■  0)  •« 
(y  .-I 
G  G  to 
o  e  c 
U  -rH  o 
a  :o  u 


a  I 

3  s 

c  o 


<  a:  5 

l-H  to 

O  <  I 


W  .  05  to 
H  o  ^  M 


(d  (U  I 

G  ti  ^ 

M-I  tJj  O'  O 

^  cn  £  ^ 
m  O  ^  § 

X  ^  ^  -u 

'  tj  • 

«  o  +j 

B  X  cd 

O  c  +>  <y 

^  c>  x: 

to  .  <M 

‘  ^  0*  -r^  O 

I  2  ^  c 

;  £  ^  <p  -2 

^  o  c  -c  W 

[  +J  ^  t»-l 

■  2  tJ  C  *'^ 

’30)0  *0 

*  O  tn  yj 

c  ^  j;  ^ 

S  S 

'  -w-t  ^  JZ  0)  c 

:  g  M  o  ® 

;  «  -s  ^  « 

<  ’-0  UJ  «  c 

>  ^  .2  o 


G  G  to 
C  E  C 
U  ♦r^  O 
c,  to  y 


a  I 

3  s: 

o  o 

G  1—1 

I 


0)  to 

<  -o  05 
pi  Qi  u- 

Ml.  cd 

^  ti.  05  G 
S  lO 
ti3  .  05  OD 
H  O  ^ 


cj  0)  •  'I 

U  U  ! 

<1-1  5J  ( 

o  ^ 

Vi  > 

0)  C  : 

to  O  ^  , 

o  •'t  5 

■  X  . 

:  y  o 

W  ^ 

•  B  2  , 

o  c  . 

1  L.  z.  O 

►  tl_i  05  ^ 

W  +J 

:  I— I 

)  X  <y 

:  ?  <“ 

!  o  °  .2  ■ 

;  ^  to 

i  V  *^  -^  ' 

:  ^  o  {:■ 

i  3  05  O 

<  o  {,_c  a 

O  ^  <y 
3  0}  OZ 

^  1  g  ^ 

3  >  5  » 

■«  O  „ 

j  x»  5 

to  to  — 

^  o  <d  5 

3  "*  CO  L, 

0  AJ  m 

0  05  -c 

3  E  'c.  ‘So 
C  E  C 
Z  u  ^  o 
-I  a  (o  u 


I  . 

s  c  ^ 
a>  c 
h-*  a 
(A  £ 
x:  p  u 
H  a; 
o 
5n 

a> 

<D  *5  c 

xi  o  <i> 


TJ  «J 

ed  c:  <i: 

^  c 

®  r-t 
4J  <1) 

c  c  cn 
o  c:  to 
m  3  o 
oj  E-i  *-• 
u 

u  -a 
c 

10 


a  c:  ^ 

nl  «  c. 

Qr  <  ^ 

o  TO 

tM  ^ 

O  CO 


CO 

o; 

c 

CO 

CO 

JC 

1 

x: 

o 

CO 

+J 

B 

CO 

ti-i 

CO 

Mi 

o 

o 

*o 

C 

u, 

Cd 

m: 

-i-t 

o 

fO 

B 

CD 

4-> 

4-[ 

■rH 

JZ 

o 

CO 

■o 

(D 

<D 

3 

B 

■r-i 

f— 1 

O 

cd 

Z 

CD 

0} 

B 

SS 

to 

SO 

<d 

o  o 

c; 

CL 

o 

CD 

B 

o 

x: 

+J 

CO 

o 

o 

+-» 

to 

u 

d-r 

-«-I 

■tH 

lO 

B 

u 

£ 

XJ 

*— 1 

-rH 

0) 

TO 

OJ 

4-> 

B 

o 

a> 

y 

c 

■4-> 

B 

cd 

o 

O  I- 
C3  TJ 


u  o  » 
o  o 

C3.  O  C 

a>  F-  ^ 

oi 

+j  ^ 
CO  cd  Q> 
•iH  Vi  B 

x:  a:  cd 

H  CL  a. 


I  tt  . 

B  c  ^ 

<1>  -rH  cT 

H  -»->  <4 
CO  E 
x:  a> 

^  pi  ^ 


4J  XI  cJ 
cd  c  < 
'll  c 


B  c  a> 

q;  B  lO 


a  xs  %« 
c  ta 

CO  X3 

u  ^  E 
<D  a; 

c,  p 
(d  2 

iD,  <  0? 

u  TO 


x: 

Vi 

<e-i 

< 

-U 

3 

O 

CO 

« 

ID 

CM 

Oi 

to 

CO 

JZ 

1 

X 

(D 

CO 

4-> 

.-1 

-P 

CkD 

B 

B. 

CO 

1 

o 

Hi: 

O 

CM 

X3 

B 

B 

CM 

b£ 

J-H 

-rH 

O 

CO 

C 

ID 

tvi 

x: 

od 

CO 

-P 

XJ 

ID 

0) 

.»r>S 

3 

B 

1-^ 

O 

cd 

Vi 

s 

0) 

<l>  c 

CO  CO  OS 

D. 

4)  "'H 
x:  -M  lO 
+J  CO  t> 

1-1  -vH 
B  Q 

W  CD  flJ 
+J  C 
<D  O  >1 
B  «d  O 

°  ^  X5 

.2  S  3 


C3  at 

o  o 

a  0}  B 

(D  B  >— 
Q£  B 

+>  »— I 

CO  od  e> 
B  c 

JB  o  cd 

H  csr  a. 


perature  Gas  Characteristics  Meet:lng  of  the  AGARD  Wind  Tunnel  and  Model  Testing  perature  Gas  Characteristics  Meeting  of  the  AGARD  Wind  Tunnel  and  Model  Testing 

Panel  (now  Fluid  Dynamics  Panel)  held  21-23  September  1959,  in  Aachen,  Germany.  Panel  (now  Fluid  Dynamics  Panel)  held  21-23  September  :959.  in  Aachen,  Germany. 


